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Abstract 



We derive the semiclassical WKB quantization condition for obtaining the en- 
ergy band edges of periodic potentials. The derivation is based on an approach 
■ which is much simpler than the usual method of interpolating with linear potentials 

in the regions of the classical turning points. The band structure of several periodic 
^ ' potentials is computed using our semiclassical quantization condition. 

§ . PACS number(s): 31.15.Gy, 31.70.Ks 

T— I \ 

§^ : 11. Introduction 

p • The study of periodic potentials is of both physical and mathematical interest. For in- 

stance, in condensed matter physics, knowledge of the existence and locations of band 
^ ! edges and band gaps of periodic potentials is very important for determining many phys- 

^ \ ical properties. Although many general mathematical properties of the eigenstates are 

^1 known unfortunately, even in one dimension, there are very few analytically solvable 

periodic potential problems. 
^ ■ One of very successful methods in numerous applications in physics and mathematics 

is the semiclassical approach. The semiclassical WKB approximation for one dimensional 
potentials with two classical turning points is discussed in most quantum mechanics text- 
books . It was originally proposed for obtaining approximate eigenvalues in the limiting 
case of large quantum numbers. It has been successfully used for many years to deter- 
mine eigenvalues and to compute barrier tunneling probabilities. The analytic properties 
of the WKB approximation have been studied in detail from a purely mathematical point 
of view, and the accuracy of the method has been tested by comparison between analytic 
and numerical results [Q, ^ There has been a special surge of interest in recent years 
due to the development of the supersymmetric WKB method f^, 0, 0. 

In this paper, we generalize the standard WKB approach to treat periodic potentials. 
We derive a new semiclassical quantization condition whose solutions give the energy band 
edges of periodic potentials. Our derivation is considerably simpler than the standard 
approach which makes use of connection formulas to match the WKB solutions in 
classically allowed regions with the solutions in classically forbidden regions. In the usual 
approach, one makes a linear approximation to the potential in the regions around the 
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classical turning points where the semiclassical first order WKB wave function diverges. 
Although the connection formulas resulting from the usual approach turn out to be simple 
enough, the derivation is quite tedious. Our simpler approach effectively amounts to 
matching the values and derivatives of the zeroth order WKB wave function which is 
non-divergent at the classical turning points 0, pJ| . 

In Sec. 2, we describe and justify our simpler approach by first re-deriving the stan- 
dard WKB quantization condition for potentials on the infinite real line with two classical 
turning points. The extension to periodic potentials of period L is given in Sec. 3, the 
main result being the quantization condition Eq. (|16[). To the best of our knowledge, 
this result and our method of derivation have not been previously discussed, even though 
the semiclassical approximation and periodic potentials have both been studied for many 
years. In order to check the accuracy of the semiclassical approach, we consider appli- 
cations to the well-studied class of Lame potentials V{x) = ma{a + 1) sn^(x,m), where 
sn(a;, m) is a Jacobi elliptic function This is a good choice, since it is one of the 



very few potentials for which the band edges are analytically known for integer values of 



13, O, A comparison of band edge energies obtained from our WKB quantization 



condition with exact results is given in Sec. 4. The limitations and successes of the WKB 
approach for periodic potentials are discussed. 



2 Simpler Derivation of the Usual WKB Condition 

In this section, we look at the standard situation of a potential on the entire real line, 
which has two classical turning points and xr given by V{x) = E for any choice of 
energy E. From now on, for simplicity, we restrict our attention to symmetric potentials 
V{x) = V{—x). For this case, xl = —xr, and it is sufficient to just look at the half 
line X > 0, since the eigenfunctions will be necessarily symmetric or antisymmetric. To 
derive the WKB quantization condition we have to connect the solution in the classically 
allowed region with the solution in the classically forbidden region. For the symmetric 
case, the zeroth order WKB approximation to the wave function for a; > is 

4'\x)=Acos[x{x)-xm (1) 
in the classically allowed region xl < x < xr and 

^W(x) = 5e-^(")+^("«) (2) 
in the classically forbidden region x > xr. We are using the notation 

1 

X(a;) = ^y p{x,E)dx, (3) 

where p{x, E) = ^2m\E — V{x)\ is the generalized momentum. Matching the wave func- 
tions ipf'\x) and ipfj{x) and their first derivatives at xr gives two equations 

Acos[x{xR)-xm = B, (4) 
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p{x,E)dx = -[n+ -y n = 0,1,2,... 



Asin[x{xR)-xm = B, (5) 
which yield tax].[x{xR) — x(0)] = 1, or 

''p(x,E)c/x = ^TT, ^7r,^7r,... (6) 

Similarly for the antisymmetric case, the zeroth order WKB approximation to the wave 
function for x > is ijjf \x) = Asm[x{x) — x(0)] the classically allowed region xl < 
X < Xji and iljfj{x) = iJe^^^^^"^-^'-^^-* in the classically forbidden region x > xr. Matching 
these wave functions and their first derivatives at xr now gives taii[x{xR) — x(0)] = — 1, 
or 

1 r^R 3 7 11 

Combining Eqs. (|]) and (|^), the quantization condition is 

% Jo 

which is the usual WKB result for a symmetric potential The derivation given above is 
evidently much simpler than the usual textbook approach for deriving connection formu- 
las. However, some comments, explanation and justification of the method used is needed 
for several points. The usual approach makes use of first order WKB wave functions 
il^^^^x) = ip''^\x) / \J p{x) ^ which diverge at the classical turning points. Although this 
divergence is understandable in the classical limit, since a classical particle has zero speed 
at the turning points, it is certainly not present in a full quantum mechanical treatment. 
Since il)'^^\x) is singular, it is necessary to resort to connection formulas and somewhat 
tricky matching of the wave function il)^^\x) and its first derivative 0, that eventually 
yields the well known WKB quantization condition Eq. (^. 

Why our simple procedure for matching il)^^\x) is justified? Clearly, the correct ap- 
proach is neither to match 'ip^^\x) nor il)'^^\x), but to keep a sufficient number of higher 
order contributions in so that the resulting wave function is non-divergent |TI|. This 
has to be the case, since there is no divergence in the full wave function. A simple 
way in which the divergence gets tamed is for the WKB wave function to have the form 
ip^^^{x) = 4'^^\x)/[^p{x) + hf{x, h)], where f{x, fi) is an analytic function of x and fi. 
It is easy to check that requiring ip^^^(x) and its derivatives to be continuous amounts 
to our procedure of matching the value and slope of at the classical turning point 

xr, which justifies our simple approach. 



3 Generalization to Periodic Potentials 

For a potential with period L, one is seeking solutions of the Schrodinger's equation 
subject to the Bloch condition 
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^P{x) = e'^^ ij{x + L) , (9) 

where k denotes the crystal momentum. The spectrum shows energy bands whose edges 
correspond to kL = 0,7r, that is the wave functions at the band edges satisfy ip{^) = 
±'?/'(x + L). For periodic potentials, the band edge energies and wave functions are often 
called eigenvalues and eigenf unctions, and we will also use this terminology. A general 
property of the eigenfunctions for a potential with period L is the oscillation theorem 
which states that the band edge wave functions arranged in order of increasing energy 
have periods L, 2L, 2L, L, L, 2L, 2L, . . .. 

For any periodic potential, it is sufficient to consider just one period of width L, say 
the interval [— f ]• In this paper, we are discussing analytic potentials with Vmin and 
Vmax as the minimum and maximum values. Further, we are only looking at symmetric 
potentials V{x) = V{—x), which necessarily makes x = a maximum or minimum. Let 
us take the origin x = to be at Vmin- The eigenfunctions will either be symmetric (S) 
or antisymmetric (A) about x = 0. Furthermore, it is easy to see that the potential is 
also symmetric about x = ^, since periodicity and symmetry about x = imply that 
V{x + f ) = V{—x + "I). Consequently, x = -I is an extremum, and the eigenfunctions 
are necessarily symmetric or antisymmetric about x = ^. Clearly, there are four types 
of eigenfunctions (S,S), (A,S), (S,A), (A, A), where the first letter denotes symmetry or 
antisymmetry about the origin, and the second letter denotes symmetry or antisymmetry 
about the point x = ^. Note that (S,S) and (A, A) wave functions have period L, whereas 
(A,S) and (S,A) wave functions have period 2L. The ground state, being nodeless, is of 
type (S,S). 

To obtain the WKB quantization condition for a periodic potential, consider any 
energy E which gives two classical turning points xl and xr in the interval [—■§,■§]• Since 
we are looking at symmetric potentials, clearly Xl = —Xr, and one only needs to look 
at the half-interval [0, ^]. Note that the region |x| < xr is classically allowed (region I) 
whereas the region xr < \x\ < ^ is classically forbidden (region II). Let us now derive 
the WKB quantization condition corresponding to the four eigenfunction types. 

For the (S,S) case, the zeroth order WKB approximation to the wave function is given 
by Eq. (|I]) in region I and 

ipfj^ix) = Bcosh[-x{x)+x{L/2)] (10) 

in region II, where x{x) is defined by Eq. (^. Matching the wave functions iljf'\x) and 
i/jfj^x) and their first derivatives at xr gives 

AcoslxixR) - xm=B cosh[-xixR) + x{L/2)], (11) 

Asm[x{xR)-xm=Bsmh[-xixR)+xiL/2)], (12) 
which, on division, give the equation 



tan 



1 /"^'^ 

p(x, E) dx 



h Jo 



tanh 



1 n 



p{x, E) dx 



(13) 
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Similarly, for the (A,S), (S,A) and (A, A) cases, we take zerotli order WKB wave 
functions with appropriate symmetry 



(AS) 



ipf\x) = Asm[x{x) - X 



4°;(x)=5cosh[-x(x) + x(V2)], 



{S,A): 



ipf\x) = Acos[x{x) - x(0)], 
= Ssinh[-x(x) + xiL/2)], 



',(0) 



{A, A): 



ipf\x) = Asm[x{x) - X 



ipf^x) = Bsinhi-xix) + x{L/2)], 
and match wave functions and their first derivatives at xr. Combining all results gives 



tan 



tan 



h Jo 



p{x, E) dx 



h Jo 



XR 



p{x, E) dx 



± tanh 



± coth 



1 

L 



p{x, E) dx 



p{x, E) dx 



Xr 



(14) 



(15) 



Taking into account the periodicity of tanx(a;), we can write a single combined equation 



h 



Xr 



p{x, E) dx 



-n ± arctan 



tanh \ t j ^ dx 



(16) 



This is our final semiclassical WKB quantization condition for symmetric periodic 

potentials of period L, the solutions to which are the WKB band edge energies E^^^. 

The quantum number n takes on non-negative integer values which keep the right hand 

side of Eq. ( |T^ ) positive. Eq. (0) is a generalization of the usual two classical turning 

points quantization condition of Eq. (§), which is found in quantum mechanics textbooks. 

r 1 — 1 

This can be readily established, since when L oo, tanh ^ Jj^^ p{x, E) dx tends to unity, 
and one gets Eq. (pT). 



4 Applications and Discussion 

In order to study the accuracy of band edges resulting from the WKB quantization con- 
dition derived in the previous section, we consider several examples. These are selected 
from the class of Lame potentials 



V{x) = ma{a + 1) sn^(a;, m). 



(17) 
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where the Jacobi elhptic function sn(x, m) has period 4K(m). The potentials have a 
period L = 2K{m). It is well known that for any integer value a = 1,2,3,..., the 
corresponding Lame potential has 2a + 1 band edges corresponding to a bound bands 



followed by a continuum. All band edge eigenfunctions are analytically known [|T^, |T5 

The Lame potential given in Eq. (^) has Vmm = and Vmax = rna{a + 1). A 
comparison of exact, analytically available band edge energies with results obtained from 
the WKB quantization condition applied to a variety of Lame potentials is shown in Table 
1. 

The quantum umber n and the symmetry of the wave functions about the points a; = 
and X = ^ are also tabulated. The results are consistent with the pattern required by 
the oscillation theorem for periodic potentials Note that for m = 1, the period of the 
Lame potentials becomes infinite and all finite bands reduce to zero width. This feature, 
which is evident for the exact energies of the potential V = 12 sn^(a;, 1) considered in Table 
1, is also reproduced in the WKB approach. Also, in Table 1, it should be noted that 
some exactly known band edges occur above Vmax- For example, choosing a = 2,m = 0.5, 
one has the the potential V{x) = 3sn^(x,0.5). It has two energy bands ranging from 
1.27 to 1.50 and from 3.0 to 4.5, with a continuum above 4.73. Clearly, the band edges 
above Vmax = 3 cannot be obtained by the semiclassical WKB method, since there are no 
classical turning points for E > Vmax- This is a limitation of the semiclassical approach. 
However, for Vmm ^ E < Vmax, it is apparent from Table 1 that results obtained from our 
WKB quantization condition are in modest agreement with exact band edge eigenenergies. 
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Table 1: Exact and WKB band edge energies. 



Potential 


^exact 


^WKB 


n 


Symmetry 


3sn2(x, 0.5) 


1.27 


1.34 





S,S 




1.50 


1.96 


1 


A,S 




3.00 


2.81 


1 


S,A 




4.50 










4.73 








6 sn^(a;, 0.5) 


2.05 


2.19 





s,s 




2.13 


2.35 


1 


A,S 




5.05 


4.95 


1 


S,A 




6.00 










6.95 










9.87 










9.95 








9.6sn^(x, 0.8) 


2.68 


2.87 





s,s 




2.68 


2.87 


1 


A,S 




7.04 


7.10 


1 


S,A 




7.20 


7.49 


2 


A,A 




9.32 


9.12 


2 


s,s 




10.52 










10.96 








12sn2(a;,l) 


3.00 


3.21 





s,s 




3.00 


3.21 


1 


A,S 




8.00 


8.14 


1 


S,A 




8.00 


8.14 


2 


A,A 




11.00 


11.07 


2 


S,S 




11.00 


11.07 


3 


A,S 




12.00 


12.07 


3 


S,A 
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